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The present paper is devoted to the analysis of cosmological solutions and its stability in the 
frame of F(R) Hofava-Lifshitz gravity. The perturbations around general spatially flat FLRW 
solutions are analyzed and it is shown that the stability of those solutions depends on the type of 
theory, i.e. on the form of the action F(R), as well as on the extra parameters contained in every 
Hofava-Lifshitz theory (due to the breaking of Lorentz invariance). The (in)stability of cosmological 
solutions can provide a constraint of the models and it may give new observational predictions. A 
natural explanation of the end of inflation and radiation/matter phases can be provided by this 
class of theories. An explicit example of F(R) gravity is studied, and the transition between the 
different epochs of the Universe history is achieved. 
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CLi" I. INTRODUCTION 

, Since observational data suggests that the Universe expansion is accelerating, a large number of models have been 
proposed to explain this phenomenon. In the frame of General Relativity (GR), the accelerating expansion can not 
be explained unless new terms or fields are considered. The main candidate for dark energy is the so-called ACDM 
model, which incorporates a cosmological constant and late-time acceleration can be achieved. Nevertheless, ACDM 
model contains several unresolved problems, as the fine tuning problem, such that other different proposals have been 
seriously considered, as for example scalar fields (quintessence, phantom,..) or modifications of GR as the so called 
F(R) gravity (see Refs. [1[ for reviews on unification of inflation and dark energy in modified gravity). At the same 
time, it is also accepted that another epoch of accelerated expansion, known as inflation, occurred during the early 
Universe. This suggests that both epochs may be unified under the same mechanism. In the frame of F(R) gravity, 
! the unification of both epochs is easily achieved and it gives a natural explanation in terms purely of gravity (see 
^ ■ Refs. drU). 

Recently, a new theory of gravity that claims to be power-counting renormalizable has been suggested in Ref. [f| . 
This new theory, already known as Hofava-Lifshitz gravity, breaks Lorentz invariance, what makes the theory to be 
renormalizable, but it produces consequently a lot of problems. However, it is conjectured that the Lorentz invariance 
is recovered in the IR limit (see Ref. [7|). Some aspects of cosmology has already been studied in the frame of this new 
• theory (see Ref. [8]). Nevertheless, as in General Relativity, Hofava-Lifshitz gravity can not explain dark energy epoch 
as well as inflation without new terms or fields, remaining such problem unresolved. An extension of the standard 
F(R) gravity to Hofava-Lifshitz theory has been performed (see Refs. [9H12|]L which seems to be also renormalizable, 
and late-time acceleration can be reproduced without cosmological constant or any other exotic field (see Refs. |^.licj). 
Even the unification of dark energy epoch and inflation can be performed in this new class of theories, and the so- 
called viable F(R) models, which avoid violations of the local gravity tests, can be easily extended to Hofava-Lifshitz 
gravity (see Ref. 

At the current paper, cosmological solutions of the type of spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) are studied in the frame of F(R) Hofava-Lifshitz gravity, and their stability is analyzed. In particular, we 
focus on the study of stability of radiation/matter dominated eras, where the Universe expands by a power law, and 
de Sitter solutions, which can well describe the accelerated expansion epochs of the Universe history. We explore 
spatially independent perturbations around these solutions, where the effects of the extra terms incorporated by the 
function F(R) are studied. Also the new parameters included in the theory, due to the breaking of the Lorentz 
invariance, could affect the cosmological solutions. The (in)stability of a solution gives very important information, 
as the possible exit from a phase of the cosmological history, constraints on the kind of action F(R) and/or future 
observational predictions. An explicit example of a F(R) action, where an unstable de Sitter solution is found, is 
studied. This model performs a successful exit from inflation, and produces an instability at the end of matter 
dominated epoch, such that a phase transition may occur. 

The paper is organized as follows: in the next section, F(R) Hofava-Lifshitz gravity is briefly introduced, and the 
cosmological equations are obtained. Sect. Ill is devoted to the analysis of general spatially flat FLRW solutions 
and their stability. The general equation for the perturbations in the linear approach is obtained. In particular, the 
class of solutions described by a scale factor that depends on a power of time (radiation/matter dominated epochs) 
are analyzed. In Appendix A, de Sitter solutions are also studied in detail. In Sect. IV an explicit example of 
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F(R) gravity is studied, where one of the so-called viable models, that unifies dark energy and inflationary epochs, is 
analyzed. Finally, some discussions and conclusions are provided in the last section. 



II. FRAMEWORK 

In this section, modified Hofava-Lifshitz F(R) gravity is briefly reviewed (see Refs. We start by writing 

a general metric in the so-called ADM decomposition in a 3 + I spacetime (for more details see Refs. [l4[ and 
references therein), 

ds 2 = -N 2 dt 2 + (dx l + N l dt)(dx j + N J dt) , (1) 

where i,j = 1,2,3, N is the so-called lapse variable, and N 1 is the shift 3-vector. In standard General Relativity, the 
Ricci scalar can be written in terms of this metric, and yields 

R = K tJ K l] -K 2 + i? (3) + 2V M KV ' v n v - n^n") , (2) 

here K = g*i Kij, Kjj is the extrinsic curvature, R^ is the spatial scalar curvature, and a unit vector perpendicular 
to a hypersurface of constant time. The extrinsic curvature Kij is defined as 

K ij = ±- T ($ ) -V^Nj-V^Ni) . (3) 
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In the original Hofava-Lifshitz model @, the lapse variable iV is taken to be just time-dependent, so that the 
projectability condition holds and by using the foliation-preserving diffeomorphisms (j^J), it can be fixed to be N = 1. 
As it is pointed out in imposing the projectability condition may cause problems with Newton's law in Hofava 
gravity. On the other hand, Hamiltonian analysis shows that the non-pro jectable f (-R)-model is inconsistent (see 
Ref. |l6(). For the non-pro jectable case, the Newton law could be restored (while keeping stability) by the "healthy" 
extension of the original Hofava gravity of Ref. [l5[ . 

The action for standard F(R) gravity can be written as 



S = J <rxyfg$)NF(R) . (4) 

Hofava-Lifshitz gravity is assumed to have different scaling properties of the space and time coordinates 

x l = bx\ t = b z t , (5) 

where z is a dynamical critical exponent that renders the theory renormalizable for z — 3 in 3 + 1 spacetime dimensions 
as it is shown in [6| (for a proposal of covariant renormalizable gravity with dynamical Lorentz symmetry breaking, 
see [l7|). GR is recovered when z = 1. The scaling properties ([5]) render the theory invariant only under the so-called 
foliation-preserving diffeomorphisms: 

fa^CCzV), ft = /(*). (6) 

It has been pointed out that, in the IR limit, full diffeomorphisms are recovered, although the mechanism for this 
transition is not physically clear. The action considered here was introduced in Ref. 

S = -^j / dtd 3 xy/g~WNF(R) , R = K {j K ij - \K 2 + i? (3) + 2 M V ' M (n"V ' v n v - n v V v n^) - L (3) ( ffl (3) ) , (7) 
2k J 

where K is the dimensionless gravitational coupling, and where, two new constants A and \x appear, which account 
for the violation of full diffcomorphism transformations. A degenerate version of the above F'(i?)-theory with /i = 
has been proposed and studied in Ref. [l2[ • Note that in the original Hofava gravity theory Q , the third term in the 
expression for R can be omitted, as it turns out to be a total derivative. The term L^(g^) is chosen to be @ 

LW(g&)=E i *G m E kl , (8) 
where the generalized De Witt metric is given by, 

Qijkt = I ^ g (3)ik g (3)jl + g(3)a g (3)jk^ _ ^ g (3)ij g (3)kl _ (gj 
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In Ref. [gj], the expression for is constructed to satisfy the "detailed balance principle" in order to restrict the 
number of free parameters of the theory, and it is defined through the variation of an action 

^ Eij = SW]gu\ (1Q) 

where the form of W[*?jm] is given in Ref. [l8[ for z = 2 and z = 3. Other forms for L^(g^) have been suggested that 
abandons the detailed balance condition but still render the theory power-counting renormalizable (see Ref. [Io|). 

We are interested in the study of cosmological solutions for the theory described by action Q. Spatially- flat FLRW 
metric is assumed 

3 

ds 2 = -N 2 dt 2 + a 2 (t)Y / (dx 1 ) 2 • (11) 
»=i 

If we also assume the projectability condition, N can be taken to be just time-dependent and, by using the foliation- 
preserving diffeomorphisms ([6]), it can be set to unity, N = 1. When we do not assume the projectability condition, 
N depends on both the time and spatial coordinates. Then, just as an assumption of the solution, N is taken to be 
unity. 

For the metric (TTT1) . the scalar R is given by 

3(1 - 3A + 6p)H 2 6(i d (H 



R = ^ n^^^jAn)- (12) 

For the action 0, and assuming the FLRW metric (fT2"j). the second FLRW equation can be obtained by varying the 

(3) 

action with respect to the spatial metric , what yields 

= F(R) - 2(1 - 3A + 3/i) (h + 3H 2 ) F'(R) - 2(1 - 3X)HRF"(R) + 2fi (r 2 F (3 \R) + RF"{R)) + n 2 Pm , (13) 



here k 2 = 16irG, p m is the pressure of a perfect fluid that fills the Universe, and N = 1. Note that this equation turns 
out the usual second FLRW equation for standard F(R) gravity (0J when A = \i = 1. If we assume the projectability 
condition, variation over N of the action ([7]) yields the following global constraint 

= J d 3 x F(R) ~ 6F'{R) | (1 - 3A + 3fi)H 2 + /iijj + 6(iHRF"(R) - n 2 p m . (14) 

Now, by using the ordinary conservation equation for the matter fluid p m + 3H(p m +p m ) = 0, and integrating Eq. (| 13[> . 
it yields 

= F(R) - 6 [(1 - 3A + 3p)H 2 + u,h] F'(R) + 6p J HRF"(R) - n 2 p m ~ , (15) 

J a 6 

where C is an integration constant, taken to be zero, according to the constraint equation (|14[) . If we do not assume 
the projectability condition, we can directly obtain (|15p . which corresponds to the first FLRW equation, by varying 
the action (J?} over N. Hence, starting from a given F(R) function, and solving Eqs. (|13p and (jl~4"]) . a cosmological 
solution can be obtained. 



III. COSMOLOGICAL SOLUTIONS AND ITS STABILITY IN F(R) GRAVITY 

In this section, we are interested to study the stability of general cosmological solutions in the frame of F(R) 
Hofava-Lifshitz gravity, with special attention to those cosmological solutions that shape the history of the Universe, 
as de Sitter or power law solutions. It is well known that in standard F(R) gravity, any cosmological solution can be 
reproduced by reconstructing the function F(R) (see Ref. [5j). As it was shown in Ref. [llj, dark energy and even the 
unification with the inflationary epoch can be reconstructed in this new frame of F(R) Hofava-Lifshitz theories. The 
stability of those solutions plays a crucial role in order to get the transition from one cosmological phase to another. 
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A. Stability of general fiat FLRW cosmological solutions 

Let us start by studying a general spatially flat FLRW metric (fTTj) . We focus specially on de Sitter and power law 
solutions of the type a(t) oc t m because dark energy and radiation/matter dominated eras are governed by this class 
of cosmological solutions respectively. The implications of the extra geometrical terms coming from F(R) could be 
determinant for the stability and transition during the epochs of Universe evolution. Firstly, we assume a general 
solution, 

H(t) = h(t) . (16) 

Then, the scalar curvature R yields, 

R h {t) = 3(1 - 3A + 6p)h 2 {t) + 6fxh(t) . (17) 

Assuming a certain F(R) that reproduces the solution (fl~6|). the FLRW equation (fl~5)) has to be fulfilled, 

= F(R h ) - 6 [(1 - 3A + 3p)h 2 + ph] F'(R h ) + 6phR h F"{R h ) - k 2 p m , (18) 

where the matter fluid is assumed to be a perfect fluid with equation of state p m = w m p m , where w m is a constant. 
By the energy conservation equation p m + 3h(l + w rn )p. m = 0, the evolution of the matter energy density can be 
expressed in terms of the solution h(t) as, 

Pmh =P^- Z{1+w - )!h{t)d \ (19) 

where po is an integration constant. We are interested to study the perturbations around the arbitrary solution h(t). 
For that purpose, let us expand the function F(R) in powers of R around (fTT)) , 

F(R) = F h + F' h {R - R h ) + -h-(R - R h ) 2 + -%-(R - R h ) 3 + 0(R - R h f , (20) 

2 o 

where the derivatives of the function F(R) are evaluated at i?/>, given in (fTTj) . Note that matter perturbations also 
contribute to the stability, inducing a mode on the perturbation. Then, we can write the perturbed solution as, 

H{t) = h(t) + 6(t) , Pm ~ p mh (l + S m {t)) . (21) 

Hence, by introducing the above quantities in the FLRW equation, the equation for the perturbation S(t) becomes 
(in the linear approximation), 

i+bS+ - H =^Mk: s - (22) 



where. 



ti l-3A + 3u 1-3A + 6m ox B v\ ^ 3) 

b= -h+ +6((l-3X + 6p)hh + ph) , 

hp p \ ' *h 

pi 

lu 2 = [1 -3A + 6a*- 2(1 - 3\ + 3p)h] ' " 



6p 2 hF£ 

o> r, x /-l + 3A-3u, -l + h ; \ h „. oA ^^/l-3X + 6p- -A F, (3) 
+ (1 - 3A + 6/i) j -h + —h + - + 6(1 - 3A + 6/i) —hh + h . (23) 

In this case, the solution for 5(t) can be split in two branches, one corresponding to the homogeneous part of the 
equation (|22[) , whose solution will depend on the background theory, i.e. on F(R) and its derivatives, and another 
one corresponding to the particular solution of eq. (|22[) . which represents the term induced by matter perturbation 
5 m . Then, the complete solution can be written as, 



5{t) = 5 homg (t) + S inh (t) 



(24) 
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We are interested in the perturbations induced by the function F(R) and its derivatives, so that we focus on the 
homogeneous solution bhomg- By a first qualitative analysis, we can see that the homogeneous part of the equation 
(|22f yields exponential or damped oscillating perturbations. The form of the perturbations depends completely on 
the form of the function F(R) and its derivatives evaluated at Rh- Note that in general, the equation (|2"3"|) has to be 
solved by numerical methods. Nevertheless, we could assume some restrictions to obtain qualitative information. Let 
us consider the cases, 

• The trivial case, given by F' h = F% = F^ = 0, makes the perturbation tends to zero, S(t) = 0, and the 
cosmological solution h(t) is stable. 

• For F' h ^ and Fjf, F^ -> 0, the term that dominates in ([22} is given by, 

uj 2 ~ [1 - 3A + 6 M - 2(1 - 3A + 3/x)fe] ^f^,, . (25) 

And the stability of the cosmological solution h(t) depends on the sign of this term, and therefore, on the model 
F(R) and the solution h(t). 

• For Fl,F£ — > but F^ 3 * 1 =/= 0, looking at (j2"3")l . the perturbation depends on the value of the last term in the 
coefficients b and ui 2 , which can be approximated to, 

b~6((l-3\ + 6u)hh+uh) , uj 2 ~ 6(1 - 3A + 6fi) | 1 - 3A + ^ hh + h ) . (26) 
V ' F h V A* / -Pfc 

The cosmological solution will be stable in the case that both coefficients (|2"6")1 are greater than zero, what yields 
a damped oscillating perturbation that decays. 

However, in general the equation (|2"21 can not be solved analytically for arbitrary solutions h(t) and actions F(R), and 
numerical analysis is required. Nevertheless, by imposing certain conditions on F(R) as above, qualitative information 
can be obtained. In order to perform a deeper analysis, some specific solutions h(t) are studied below, as well as an 
explicit example of F(R). 

B. Stability of radiation/matter eras: Power law solutions 

In this section, an important class of cosmological solutions is considered, the power law solutions, which are 
described by the Hubble parameter, 

Tn 

H{t) = — -> a(t) oc t m . (27) 

In the context of General Relativity, this class of solutions are generated by a perfect fluid with equation of state 
parameter w = — 1 + and the matter /radiation dominated epochs are approximately described by (f27|) . Also 
phantom epochs can be described by this class of solutions when m < . Let us study the stability for the Hubble 
parameter (|27[) . and how the inclusion of extra terms in the action and the new parameters (A, fx) may affect the 
stability of the solution (|27|) . As in the above section, the perturbation equation (|22|) can not be solved analytically 
in general, although under some restrictions we can obtain important qualitative information about the stability of 
the solution. Then, by assuming an F(R) that approximately does not deviate from Hilbert-Einstein action during 

radiation/matter dominated epochs, the second and third derivatives can be neglected F^F^ ~ (as they must 
become important only during dark energy epoch and/or inflation). In such a case, the coefficient in front of S(t) in 
the eq. (|22|) is approximated as, 

u 2 ~ [1 - 3A + 6^ - 2(1 - 3A + 3/*)M«)] ^fipT, ■ ( 28 ) 

Then, the value of the frequency uj 2 depends on the time, such that the stability may change along the phase. For 
small values of t, the frequency takes the form lu 2 ~ —2(1 — 3A + 3//) 6/ tv"f h , and assuming A ~ the perturbations 

h 

F' 

will grow exponentially when -^k > 0, and the solution becomes unstable. While for large t, the frequency can be 



occur. 



approximated as uj 2 ~ (1 — 3A + 6/x) h^t)F" anc ^ ^ ne instability will be large if < 0, and a phase transition may 
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IV. EXAMPLE OF A VIABLE F(R) MODEL 



Let us consider an explicit model of F(R) gravity in order to apply the analysis about the stability performed 
above. We are interested to study the stability of radiation/matter dominated eras as well as de Sitter solutions for 
an explicit F(R). Here we consider a model proposed in Ref. @, and studied in Ref. Q in the context of standard 
gravity and generalized to Hofava-Lifshitz gravity in Ref. [ll|. The action is defined as, 

where (x, a, j3, 7) is a set of constant parameters of the theory. In standard gravity, this model can reproduce well 
late-time acceleration with no need of a cosmological constant or any kind of exotic field, as well as also inflation, 
such that the unification of both epochs of the Universe history under the same mechanism can be performed (see 
Ref. [2]). For simplicity, we assume n = 2 in (129!) for our analysis. The radiation/matter dominated epochs, which 
can be described by the class of solutions given in (|2T[) , could suffer a phase transition to the era of dark energy due 
to the instabilities caused by the second term of the action (|2"5|) . Then, we are interested to study the possible effects 
produced by the presence of these extra geometric terms during the cosmological evolution. By assuming the solution 
(|2T|) . and following the steps described in the above section, the stability is affected by the derivatives of the function 
(|29|) evaluated in h(t) = m/t. We are interesting in large times, when the end of matter dominated epoch has to 
occur. At that moment the derivatives of F(R) can be approximated as, 

K^X, Fh^-W, F^^O. (30) 

Here for simplicity, we have assumed < j3 « 1. By means of the analysis performed in the previous section, we can 
conclude that the linear perturbation S(t) grows exponentially, and the radiation/matter dominated phase becomes 
unstable for large times, what may produce the transition to another different phase. Then, the F(R) function (|2T))) 
can explain perfectly the end of matter dominated epoch with no need of the presence of a cosmological constant. 



Let us now study the stability of de Sitter solutions (for more details on de Sitter solutions and its stability, see 
Appendix A). It is known that the model (|29p may contain several de Sitter solutions (see Ref. [ll[ and Q), solutions 
of the first FLRW equation that now turns out an algebraic equation given by, 



R n (aR n ~B) 6# 2 (-I + 3A-3m) 1 + nay Iff + i?o ( 2 7#o - nff) + R z n - L {y 2 R + 2na) 

7~> 1 \ ' / 1 

Kq - 



l+yR% (1 + 7^) S 



= 0. (31) 



This equation has to be solved numerically, even for the simple case studied here, n = 2. Nevertheless, one of the de 
Sitter points from the model (1291) is defined by a minimum of the second term in the action (|29[) . By assuming the 
constraint on the parameters py/a 3> 1, the minimum that represents a de Sitter point is given by, 

R ~ (JL) ' , F'(R ) = X , F(R ) = R Q -2A, where A ~ L . (32) 

Then, by evaluating the derivatives of (|29|) around i?o and by the equation (|A6j) . the perturbation S(t) can be 

calculated. Note that the stability condition for de Sitter solution, given by j& > 12-ffp in Appendix A, is not 

satisfied for this case as Fq >> Fq, such that the de Sitter point (|3"2"j) is unstable. By resolving eq. (| A6|) . the 
perturbation is given by exponential functions, 

5{t)=C 1 ^ t + G 2 e a - t , with a± = g o(l-3A + 3/,) 

Hence, the model (|29l) is unstable around this de Sitter point (|3"2")) , what predicts the exit from an accelerated phase 
in the near future, providing a natural explanation about the end of the inflationary epoch, or a future prediction 
about the end of dark energy era. However, the theory described by (|29p may contain more de Sitter points, given 
by the roots of equation (I31[) . which may be stable. Then, a deeper analysis has to be performed to study the entire 
Universe evolution for this model of F(R) gravity. 
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V. DISCUSSIONS 



At the present paper, we have analyzed spatially flat FLRW cosmology for nonlinear Hofava-Lifshitz gravity. 
Basically we have extended standard F(R) gravity to Hofava-Lifshitz theory, which reduces to the first one in the 
IR limit (where we assume that the parameters (A,/i) are reduced to unity). The stability of this general class of 
solutions has been studied and it is shown that it depends mainly on the choice of the function F(R) and in part on 
the values of the parameters (A,/i). For large times, when the scalar curvature is very small, the main effect of the 
perturbation on a cosmological solution is caused by the value of the derivatives of F(R). It is shown that in general, 
the perturbation equation can not be solved analytically, even in the linear approach. Nevertheless, under some 
restrictions, important information is obtained, and the (in)stability of the different phases of the Universe history 
can be studied. For specific values of the derivatives of F(R), a given solution can becomes (un)stable, which means 
a major constraint on models. By analyzing an explicit example in Sect. IV, where an F(R) function of the class of 
viable models is considered, we have found that this kind of theories can well explain the end of matter dominated 
epoch, and reproduces late-time acceleration. We have shown that for this specific example, there is a de Sitter point 
that becomes unstable, what predicts the end of this de Sitter epoch, providing a natural explanation of the end of 
inflationary era. However, as this model (and in general every F(R) model) may contain several de Sitter solutions, 
where some of them can be stable, a further analysis of the phase space has to be performed to connect the different 
regions of the Universe history. 

Hence, the analysis made here provides a general approach for the study of spatially flat FLRW solutions in the frame 
of higher order Hofava-Lifshitz gravities, which can restrict the class of functions F(R) allowed by the observations, 
and it gives a natural explanation of the end of inflation and matter dominated epoch, shaping the Universe history 
in a natural way. 
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Appendix A: de Sitter solutions in F(R) gravity 

Let us consider one of the simplest but most important solutions in cosmology, de Sitter (dS) solution. As dark 
energy and inflation can be shaped (in its simplest form) by a dS solution, its stability becomes very important, 
specially in the case of inflation, where a successful exit is needed to end the accelerated phase occurred during the 
early Universe. In general, standard F(R.) gravity contains several de Sitter points, which represent critical points 
(see ||). The analysis can be extended to F(R) Hofava-Lifshitz gravity, where the de Sitter solution H(t) = Ho has 
to satisfy the first equation FLRW equation (TT5]) , 

= F(Rq) - 6ff 2 (l - 3A + 3[i)F'(R ) , (Al) 

where we have taken C = and assumed absence of any kind of matter. For this case, the scalar R is given by, 

7? = 3(l-3A + 6Ai) J ffo ■ (A2) 

Then, the positive roots of equation (|A1[) are de Sitter points allowed by a particular choice of a F(R) function. By 
assuming a de Sitter solution, we expand F(R) as a series of powers of the scalar R around i?o, 

TP" p( 3 ) 

F(R) = F a + F>(R- R ) Ro) 2 + £ ^~(R~ ^>) 3 + O^) . (A3) 

Here, the primes denote derivatives respect R, while the subscript means that the function F(R) and its derivatives 
are evaluated at R = Rq. Then, by perturbing the solution, the Hubble parameter can be writing as, 

H(t) = H + 5{t) . (A4) 

Using (|A3j) . and the perturbed solution (|A4j) in the first FLRW equation (fl"5|) . the equation for the perturbation 
yields, 

= ^F - 3Hg(l - 3A + 3/i) - 3H Q [((1 - 3A)-F^ + 6FgH$(-l + 3A - 6/i)(-l + 3A - 3/z)) S(t) 
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-6F£nH (-l + 3A - 3/x)j(i) - 12F^ r 5[t) 



(A5) 



Here we have restricted the analysis to the linear approximation on 8 and its derivatives. Note that the first two 
terms in the equation (|A5j) can be removed because of equation (|Alj) , which is assumed to be satisfied, and equation 
151) can be rewritten in a more convenient form as, 



(3A - 1)§ - 6ff 2 (l - 3A + 6a«)(1 - 3A + 3 m) 



S(t) = 



(A6) 



Then, the perturbations on a dS solution will depend completely on the model, specifically on the derivatives of F(R), 
as well as on the parameters (A,/z). The instability becomes large if the term in front of 5(t) (the frequency) in 
the equation (|A6[) becomes negative and the perturbation grows exponentially, while if we have a positive frequency, 
the perturbation behaves as a damped harmonic oscillator. During dark energy epoch, when the scalar curvature is 
very small, the IR limit of the theory can be assumed, where A = fx ~ 1, and the frequency depends completely on 

F' F' 9 

the value of j§r. In order to avoid large instabilities during the dark energy phase, the condition j$r > 12-ffg has 
to be fulfilled. Nevertheless, when the scalar curvature is large, the IR limit is not a convenient approach, and the 
perturbation depends also on the values of (A,//). If we assume a very small Fq , the frequency in the equation (|A6j) 
dominates compared to the other terms, and by assuming A > 1/3, the stability of the solution will depend on the 

f' 

sign of -F&, being stable when it is positive. 
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